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The propagation of an electromagnetic wave in a one-dimensional fractal object, the Cantor 
set, is studied. The transfer matrix of the wave amplitude is formulated and its renormalization 
transformation is analyzed. The focus is on resonant states in the Cantor set. In Cantor sets 
of higher generations, some of the resonant states closely approach the real axis of the wave 
number, leaving between them a wide region free of resonant states. As a result, wide regions 
of nearly total reflection appear with sharp peaks of the transmission coefficient beside them. 
It is also revealed that the electromagnetic wave is strongly enhanced and localized in the 
cavity of the Cantor set near the resonant frequency. The enhancement factor of the wave 
amplitude at the resonant frequency is approximately 6/ |7j r |, where n T is the imaginary part of 
the corresponding resonant eigenvalue. For example, a resonant state of the lifetime r r = 4.3ms 
and of the enhancement factor M = 7.8 x 10 7 is found at the resonant frequency ui r — 367GHz 
for the Cantor set of the fourth generation of length L — 10cm made of a medium of the 
dielectric constant e = 10. 
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1. Introduction 

We study in the present paper the propagation of an 
electromagnetic wave in a fractal object. Takeda et al. re- 
cently reported an interesting experiment 1 in which they 
put an electromagnetic microwave into Menger sponges 
made of a dielectric substance (epoxy resin). The Menger 
sponge 2 is a fractal object of the similarity dimension 
Iog20/log3 ~ 2.7, embedded in three spatial dimen- 
sions. They set the fractal medium in a square tube, put 
an incident wave from one end of the tube, and mea- 
sured the reflection wave on the same end as well as 
the transmission wave on the other end. They observed 
at lu = 12.8GHz narrow dips both in the transmission 
amplitude and the reflection amplitude and claimed a 
strong localization of light at the frequency. We show in 
the present paper that: 

(1) the strong localization of light indeed occurs in a 
fractal object at various frequencies; 

(2) resonant states with long lifetimes are responsible 
for the strong localization. 

We thereby suggest that, in their experiment, the mi- 
crowave of frequency 12.8GHz was diverted to all direc- 
tions because of resonant scattering and hence the am- 
plitudes at the both ends of the tube were lost. 

More specifically, we theoretically calculate the wave 
propagation in the Cantor set. The Cantor set 2 (Fig. 1) 
is another fractal object of the similarity dimension 
log 2/ log 3 ~ 0.63, embedded in one spatial dimension. 
In order to realize the Cantor set experimentally, one can 
use layers of sheets of dielectric media and put the light 
normal to the sheets. 

The transmission and reflection coefficients of electro- 
magnetic waves in the Cantor set was first computed by 
Sun and Jaggard in 1991 with scattering matrices, 3 later 
by Bertolotti et al. in 1996 with transfer matrices 4 and 

*E-mail address: hatano@iis.u-tokyo.ac.jp 



Generation 
zeroth (v=0) 
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fourth (v=4) 



Fig. 1. The recursive construction of the Cantor set. We first pre- 
pare a line segment of length L, which we refer to as the zeroth 
generation (y = 0) of the Cantor set. The first generation [y = 1) 
of the Cantor set is constructed by dividing the zeroth generation 
into three segments and removing the central one. The second 
generation (v = 2) is constructed by dividing each segment of 
the first generation into three segments and removing the central 
one of each. The Cantor sets of higher generations are thus con- 
structed recursively. In the strict sense of the term, the Cantor 
set as a fractal object is the one obtained in the limit v — * oo; 
the objects of finite v should be called prefractals. We here use 
the term rather loosely and refer to the object even of a finite v 
as a Cantor set of the !/th generation. 



recently by Yamanaka and Kohmoto again with trans- 
fer matrices. 5 The novel point of the present paper is 
to attribute peaks and dips of the transmission and re- 
flection coefficients to resonance poles in the complex 
wave-number plane; to our knowledge, resonant states of 
electromagnetic waves in the Cantor set have never been 
analyzed. Ohtaka 7,8 recently derived from the scattering 
matrix of a finite system, the density of states, which is 
closely related to the resonance poles. It may be inter- 
esting to apply his formula to a Cantor set. 

We first in §2 formulate the transfer-matrix method of 
computing wave amplitudes. The transfer matrix for a 
Cantor set of the 2/th generation, T^(£), is defined by 
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the transfer matrix for a Cantor set of the (y — l)th gen- 
eration, T^-^iC), where C denotes the wave number of 
the incident light, made dimensionless by the system size 
L. Thus we obtain the transfer matrix for a Cantor set 
of an arbitrary generation recursively. Incidentally, it is 
common 9 in analyzing photonic crystals to use the scat- 
tering matrix instead of the transfer matrix; the transfer 
matrix, when we multiply it by itself many times, can be 
numerically divergent. In the present case, however, we 
multiply transfer matrices hierarchically and hence the 
numerical instability does not occur. 

Once the transfer matrix is given, we compute the 
transmission coefficient T and the reflection coefficient 
R from matrix elements of the transfer matrix as 

no 



R(0 



(t(%)) : 



T x 



(1.1) 



(1.2) 



We also define the resonance poles Cr = 6 + ' lf ]i as the 
zeros of the lower right element of the transfer matrix: 



(T<">(Cr)) 



= 0. 



(1.3) 



In §3, we demonstrate by numerical calculation that 
the reflection coefficient R and the transmission coeffi- 
cient T oscillate as we vary the frequency of the incident 
wave (Fig. 5 below). The oscillation is very rapid par- 
ticularly for media of large dielectric constants and for 
Cantor sets of higher generations. 3-5 

We also demonstrate that the resonance poles Cr = 
£ r + i?7 r appear intermittently along the real axis of the 
(dimensionless) wave number £; these resonant states are 
responsible for the rapid oscillation of the reflection and 
transmission coefficients. Each peak of the transmission 
coefficient (and hence each dip of the reflection coeffi- 
cient) corresponds to a resonance pole nearby on the 
complex wave-number plane. 

The merit of considering the resonant state is the use 
of the imaginary part of the resonant eigenvalue, rj T . The 
inverse of the imaginary part gives the resonant lifetime 
in the form 



L/ci 



(1.4) 



where L is the entire length of the Cantor set and c\ 
is the speed of light outside the Cantor set. If we put 
an incident wave packet of nearly the resonant frequency 
cj r = £ r x Ci/L, the wave packet stays inside the cavity 
of the Cantor set for as long as the resonant lifetime r r . 

We also reveal that the inverse of the imaginary part 
of the resonant eigenvalue has another indication for 
the standing wave. The standing wave resonates and is 
strongly enhanced inside the cavity of the Cantor set. 
The enhancement factor M (the magnitude of the stand- 
ing wave inside the cavity divided by the magnitude of 
the incident wave) is shown to be 

Some of the resonant states have very long lifetimes, 
again particularly for media of large dielectric constants 



1 



H H 



Fig. 2. A wave with the wave number 3"m7r/2L (m is a positive 
integer) resonates in every cavity of the Cantor set. 



and for Cantor sets of higher generations. For example, 
a Cantor set of the fourth generation v = 4 with the 
dielectric constant e = 10 and of length L = 10cm has 
a resonant state of lifetime r r = 4.3ms at the resonant 
frequency uj r — 367GHz (Table I below). For the wave 
number near such resonant states, the wave amplitude is 
strongly enhanced (or localized) inside the cavity of the 
Cantor set (Fig. 7 below). The enhancement factor M 
in the above example is no less than 10 7 , a surprisingly 
large number. 

This enhancement is a consequence of repeated re- 
flection and interference inside the cavity of the Can- 
tor set. Roughly speaking, a wave with the wave number 
3 u mir/2L (m is a positive integer) forms a standing wave 
in every cavity of the Cantor set and hence has strong 
positive interference; see Fig. 2. 

In the rest of the paper, we investigate why the reflec- 
tion and transmission coefficients oscillate rapidly and 
why some of the resonant states have very long lifetimes. 
For the purpose, we introduce in §4 a renormalization 
transformation of the transfer matrix. We show that the 
transfer matrix for the Cantor set of any generation is 
expressed in the form 



TM(C) =e iC(s< 



(1.6) 



where s^\() and r^ u \() are real functions and a z and 
u x are the Pauli matrices. With this expression (1.6), we 
renormalize the fractal structure into a uniform medium, 
within which transmission and reflection occur 
continually. 

We then analyze in §5 the renormalization transforma- 
tion 



) T .("-l) 



(«(">, r(">) 



(1.7) 



perturbatively for a small dielectric constant e ~ 1. We 
observe at the dimensionless wave number £ — 3 v rmr/2 
(m is a positive integer), a rapid increase of the reflec- 
tion part r("\ and consequently a rapid increase of the 
reflection coefficient R in the form 



R~ 2 2v \ 2 , 



(l.S 



where A = (e — l)/2 is the perturbation parameter. At 
other wave numbers, the reflection part decreases as 
~ (2/3) y and the system is renormalized to a trivial fixed 
point 



0. 



(1.9) 



We have the perfect transmission, T = 1 and R = 0, at 
this trivial fixed point. 

We finally analyze how the resonance poles ( r move in 
the complex wave- number plane as we renormalize many 
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X 



Resonance poles 



Fig. 3. For v = 0, the resonance poles are regularly aligned on 
the lower half plane of f . As we progress to higher generations of 
the Cantor set, some of the poles move toward trivial fixed points 
(indicated by crosses) on the real axis, opening at the same time 
a wide gap between them. 



(a) 



(b) 



Fig. 4. (a) A layer of the medium 2 (the gray area) sandwiched 
by layers of the medium 1. (b) Two layers of the medium 2 (the 
gray areas) sandwiched by layers of the medium 1. 



times, or as we progress to higher generations of the Can- 
tor set. We conclude in two opposite limits (the weakly 
renormalized limit |s^| 3> Ir^l in §5 and the strongly 
renormalized limit \s^ v '\ <C \r^ v '\ in §7) that, as the 
renormalization progresses, some of the resonant states 
approach trivial fixed points (1.9) (examined in §8) on 
the real axis of the wave number, leaving between them 
a wide region free of resonant states (Fig. 3). As a result, 
a wide region of nearly total reflection (R ~ 1) appears 
with sharp peaks of the transmission coefficient (T ~ 1) 
beside them. This is the main point of the present paper; 
the fractal structure of the medium is responsible for the 
behavior shown in Fig. 3. 

Incidentally, we heavily use M. Suzuki's quantum anal- 
ysis in the perturbation analyses in §§5-7. We give details 
of the calculation in Appendix A. 

2. Transfer matrix 

In this section, we formulate the transfer matrix of 
wave amplitudes in Cantor sets. We first construct the 
transfer matrix for a dielectric layer, and then for Cantor 
sets of higher generations recursively. 

Let us begin with solving the situation in Fig. 4 (a). 
Using the Lorentz gauge, we write the Maxwell equation 
in the form 

d 2 - 

AA = e ifH —A, (2.1) 



dt 2 



where Si and \ii denote the dielectric constant and the 
magnetic permeability of the medium i(— 1,2), respec- 



tively, with the light speed 



(2.2) 



Since the two transverse components of the vector po- 
tential A are decoupled, we hereafter consider only one 
component Ae lujt . Thus we have 

d 2 

—^A = -e lf i t Lu 2 A. (2.3) 
The solution of Eq. (2.3) is given by 



A(x) 



with kj 



ae lklX - 


f /3e 


-ikix 


pe ik 2 x _| 


- qc~ 


i 




h 5e~ 


-\k\x 



for x < X, 

for X < x < Y, 

for x > Y, 



fejJTiUj — uj/ci. Note that 



C2 



(2.4) 



(2.5) 



where n stands for the refraction constant. The boundary 
conditions at x — X are given by 



A(X - 0) 



1 



A'(X - 0) 



(2.6) 

A'(X + Q) (2.7) 



A(X + 0), 
1 

Mi ' M2 
and those at x = Y are given likewise. We thus have 



ae 



qe 



a—e lklX 
Mi 



Mi 



M2 



or, in a matrix form, 



V2 V ie 



-ie 



M2 



(2.8) 



(2.9) 



_L ( gi 

V2\ «' 



k 2 X 
ik 2 X 



-ie 



-ik 2 X 
-ik 2 X 



where 



x — 

fa M2 



g2Ml 

£iM2 



(2.10) 



(2.11) 



In reality, the magnetic permeability is almost constant 
for most media, 



/ii ~ /i2, and hence n 



where 



£ = — 
~ El 



is the relative dielectric constant. 
Applying the matrix 



= u~ 1 



(2.12) 
(2.13) 

(2.14) 
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lc -ik 2 x ) { q I ( 2 - 15 ) 



where the angle <j> is defined by 



e^ = n'. 



(2.16) 



(The constant <p is m fact a purely imaginary number.) 
We now express Eq. (2.15) in terms of the Pauli matrices 

as 

QlklXa * ( 5 ) = V^e-^/V^' ( J ) ■ (2-17) 
Likewise, we have 

eifel1 ^ ( ij ) = Vrie-^/W^ ( £ ) (2.18) 

for the boundary conditions at £ = Y\ Solving Eq. (2.17) 
with respect to (p,iq) and substituting the solution for 
(p,iq) in Eq. (2.18), we have 



ifciYcr, ( 7 

i<5 



i(j>a v /2 ik 2 (Y-X)a^ i<pa v /2 



e "f"»/'e 



x e 



i/3 



(2.19) 



ginfci (Y — X)(cr z cos (^+(T K sin 0) 



x e ifclX ^ ( J j • ( 2 ' 2 °) 
Hence we define the transfer matrix for the medium 2 as 

72(C) = e _i ^" Ts/ / 2 e i ^ nCTz C i ^" :r! '/ 2 = e i C n ( er ^ cos <t>+^ sin <t>) 

(2.21) 

The coefficients of the Pauli matrices in the exponent are 

"» s *=KI + ?t)4< £+1 )' < 2 ' 22 » 

The situation in Fig. 4 (b) is then described by 
K \ --T 2 (k 1 (Y / - X'))e iklX '°* f 7 



IP 



= T 2 (fc 1 (y / -X'))e 1 — e 



15 



= r 2 (fci(y' - x'))c ifci(x '- y) ^r 2 (fc 1 (r - x)) 



xe 



i/3 



(2.24) 



Hence wc define the transfer matrix for the medium 1 as 
71(0=0^. (2.25) 



We are now in the position of writing down the transfer 
matrix for the Cantor set. The transfer matrix for the 
Cantor set of the first generation is given by 

r (1) (C) = T 2 (0ri(0T 2 (^0, (2.26) 



where 



C = hL 



(2.27) 



is a dimcnsionless wave number with L being the length 
of the Cantor set. We thereby generate the transfer ma- 
trix for Cantor sets of higher generations recursively as 



->("+!) 



(0 = T {v) 



(2.28) 



where v denotes the generation of the Cantor set. We 
also define 



T(°)(C) = T 2 (C), 



(2.29) 



which makes Eq. (2.26) a part of the recursion rela- 
tion (2.28). 

Finally, we mention two symmetries of the transfer ma- 
trix. First, because the Pauli matrices cr z and a x are both 
symmetric, the transfer matrices for the media 1 and 2, 
Eqs. (2.21) and (2.25), are also symmetric. The transfer 
matrices for the Cantor sets of all generations are there- 
fore symmetric as well: 

t T H =T H_ (2.30) 

Next, we note the fact that det e arT " =1 (pi = x, y, z) for 
an arbitrary coefficient a. Because the transfer matrix 
is written as a product of exponential operators of 
the Pauli matrices, wc conclude that 



detT^ = 1. 



(2.31) 



3. Transmission coefficient, reflection coefficient 
and resonance poles 

To summarize the derivation in the previous section, 
we have for a Cantor set of the ^th generation, the rela- 
tion 

^Lout 



^4-Rout 
^Rin 



T (v \C) 



(3.1) 



where Aun and AL OUt are the amplitudes of the incoming 
and outgoing waves on the left of the Cantor set, respec- 
tively, and A R ; n and A Rout are those on the right, respec- 
tively. (We here included phase factors such as e lklXcr * in 
the wave amplitudes.) In the present section, we give the 
transmission and reflection coefficients as well as the res- 
onance condition in terms of the elements of the transfer 
matrix T'"'. 

Rearranging the transfer matrix in the form of a scat- 
tering matrix, we have 

^Lout 
^Rout 



J- 22 

1 
^22 



-T21 

2iiT 22 — Ti 2 T 2 i 

-T12 
1 



1 

Ti 2 



1 

T12 

A L in 
^Rin 



^Lin 
^Rin 



(3.2) 
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Fig. 5. The reflection coefficient R (the soiid curves) and the 
transmission coefficient T (the broken curves) are piotted for 
< v < 4 (from the top down to the bottom) and for (a) 
£ = 1.01, (c) e = 1.5 and (e) e = 10, with the resonance poles 
C = Cr = Cr + ir/r shown for (b) e = 1.01, (d) e = 1.5 and (f) 
e = 10. On the panel (a), refer to the left axis for the reflection 
coefficient and to the right axis for the transmission coefficient. 
On the panels (b), (d) and (e), the dots indicate the resonance 
poles. Note that the horizontal axis of each panel is scaled by 
the factor 3". 



where we abbreviated (TM(£)) to for simplicity. 
(We used the symmetries (2.30) and (2.31) in moving 
from the second line to the third line of Eq. (3.2).) 
We thereby obtain the transmission and reflection co- 
efficients as 

2 



1 



I 22 



and 



R 



T12 



(3.3) 



Here, the flux conservation T 
identity 



T22 

R = 1 is followed by the 



-12 1 



'-22 1 



for C e 



(3.4) 



We computed the transmission and reflection coefficients 
according to Eq. (3.3) and plotted them in Fig. 5 (a), (c) 



Fig. 5. Continued. 



and (e). Notable is the wide regions of nearly total re- 
flection (i? <~ 1) and sharp peaks of the transmission 
coefficient in the panel (e). 

We can relate the peaks and the dips of the transmis- 
sion and reflection coefficients to resonance poles. There 
are various ways of defining the resonant state. One way 
of defining the resonant state is to seek a pole of the 
transmission and reflection coefficients (3.3). The reso- 
nance condition is thus given by 



T 22 = 0. 



(3.5) 



Another way of defining the resonant state is to put the 
incoming wave amplitudes to zero; 10-12 in other words, 
the resonant state as a solution of the stationary wave 
equation is a function with outgoing waves only. (This is 
sometimes called the Siegert condition.) We obtain the 
solution by putting Ayh u — Au n = in Eq. (3.1). This 
also yields the condition (3.5). 

The condition has resonance solutions 



C = Cr = 6 + i?7r 



(3.6) 
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Fig. 5. Continued. 



Fig. 5. Continued. 



in the complex £ plane. (Note, however, that the imagi- 
nary part of the resonant eigenvalue is in fact negative: 
r\ x < 0. 11,12 A solution with a positive imaginary part 
would be a function with incoming waves only.) The res- 
onant frequency iv r is given by 



c i c 



(3.7) 



(3.8) 



and the resonant lifetime is given by 

L 1 
ci W 

We plot numerically obtained resonance solutions in 
Fig. 5 (b), (d) and (f); each dot indicates the divergence 
of 1/ |T*22 1 j or a solution of the resonance condition (3.5). 

For s — 10 and v = 4, we found resonance solutions 
with quite long lifetimes. We list in Table I the solutions 
whose imaginary parts are less than 1.0 x 10~ 5 . (We ob- 
tained the resonance solutions as follows: we first wrote 
down the resonance condition (3.5) with Mathematica 
explicitly; we then solved the equation on the complex 
wave-number plane numerically by the Newton-Lapson 



Table I. Some of the resonance poles f r = £ r + irj^ for e = 10 and 
v = 4 with particularly long lifetimes. We put L = 10cm and 
ci = 2.9979 X 10 m for calculation of the resonant frequency lu t 
and the resonant lifetime r r . 
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-2.34999 > 
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(b) 
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-7.68867 > 


< 10" 
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0.367 


4.34 


(c) 


130.988 


-5.91115 x 10" 
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0.393 
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189.149 


-8.38728 x 10" 


-« 


0.567 


0.0398 
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292.234 


-6.21849 x 10" 
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0.876 


0.0536 
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-1.43720 x 10" 


i> 


1.08 


0.232 
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-5.26455 > 
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00 
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ti 
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0.120 
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-ti 
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0.0485 
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0.0440 
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-7.07779 x 10" 


-7 


2.74 
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Fig. 5. Continued. 



Fig. 5. Continued. 



method. We can easily choose the initial guess for the 
Newton-Lapson method, since the distance between the 
real parts of the neighboring resonance poles is roughly 
equal for all poles.) 

At the resonance solutions on the complex ( plane, 
the transmission and reflection coefficients (3.3) diverge. 
The tail of the divergence forms structures on the real 
£ axis as is exemplified in Fig. 5. On the real axis near 
the resonance solutions, the transmission coefficient has 
a peak: T <~ 1; sec, for example, the resonance pole at 
(£ r ,r? r ) ~ (57r/16, -0.2) for v = in Fig. 5 (f) and the 
corresponding peak of the transmission coefficient T for 
v = in Fig. 5 (e). Because of the conservation law, the 
reflection coefficient then must have a dip: R <~ 0. 

Finally, we show the spatial distribution of the field 
A(x). We can compute the field amplitude inside the 
Cantor set by expressing the transfer matrix with 
the elementary transfer matrices T\ and T 2 explicitly. We 
apply the elementary transfer matrices one by one to the 
amplitude vector and obtain the field amplitudes inside 
the Cantor set. 



Figure 6 shows the field distribution for e = 10 and, in 
particular, for the cases that the reflection coefficient is 
close to unity. In this figure, j4 r ; g ht and A\ e ft denote the 
amplitudes of the right-going component and the left- 
going component, respectively; that is, the total field is 
given by 

(3.9) 



A{x) - A right e ik > x + A lc{t e- ik * x 



in the medium i(= 1,2). We see that the incoming wave 
from the left is bounced almost completely and barely 
transmits to the right. 

On the other hand, Fig. 7 shows the field distribution 
in the cases close to the resonance poles given in Ta- 
ble I (a) and (b). We see prominent localization of the 
field inside the Cantor set. 

Let us show that the enhancement of the field inside 
the Cantor set is of the order of the inverse of the imag- 
inary part of the resonance pole. Suppose that the inci- 
dent wave ^Lin from the left of the Cantor set of the ^ th 
generation causes a reflection wave ^Lout and a trans- 
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Fig. 6. Semi-logarithmic plots of the field distribution A(x) for 
e = 10 and < v < 4 (from the top down to the bottom) at 
£ = fciL = 3"7r/2. (The condition ( = 3"m.7r/2 shown in Fig. 2 
is on resonance only for e ~ 1. It is off resonance for e = 10.) 
Each solid line indicates the magnitude of the right-going wave, 
bright | i each broken line indicates the magnitude of the left- 
going wave, |Ai c f t | 2 , and each gray line indicates the magnitude 
of the total standing wave, |A(a;)| 2 = \ A right c lkiX + Ai e f t e _lfci:l! | 2 
in the medium i(= 1,2). 




Fig. 7. Semi-logarithmic plots of the field distribution A(x) for 
e = 10 and v = 4 at (a) ( = kiL = 47.29458732802431 and 

(b) C = kiL = 122.4274149967578. Each solid line indicates 

l 2 

the magnitude of the right-going wave, A^g^t , each broken 
line indicates the magnitude of the left-going wave, |Ai c f t | 2 , and 
each gray line indicates the magnitude of the total standing wave, 
|A right e ifc ^+A lr ~e- ifc ^l 2 



\A(xr 



iloftC 



in the medium i(= 1, 2). 



mission wave Aptout 

Aptout 





. Then we have the relation (3.1), or 
= T M(C)f ^ Lin V (3.10) 

\ ^Lout / 



Near the resonance but on the real ( axis, the transmis- 
sion coefficient is close to unity, which means | Agouti — 
\Aun\ and ALout — 0. We are now interested in the right- 
going field A r ight and the left-going field Ai e f t inside the 
cavity L/3 < x < 2L/3. This is related to A^ out as 



^Rout 





T (u-1) I C 

3 



bright 
Aloft 



(3.11) 



Let us for the moment use the following symbols for the 
sake of simplicity: 



tM(c) 

t("-i) ( i 
3 





Tl 2 


T 2 i 


T22 


in 


tl2 


*21 


*22 



(3.12) 
(3.13) 



Equation (3.11) takes the form of scattering of an in- 
coming wave A r ight due to the part 2T/3 < x < L of the 
medium 2. Hence the transmission coefficient due to the 
part 2L/3 < x < L gives the enhancement of the wave 
amplitude as follows: 

2 



I A 



Rout I 



1 



^22 



iright I 



(3.14) 
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or 



Aright | 



14 



right I 



iRout | 



1*22 r 



(3.15) 



Because of the spatial symmetry at the resonance, we 
also have | -Aright | — |4eft|- This is indeed confirmed by 
considering a relation from Eq. (3.11), 



lAwtl 2 = 



'12 



*22 



iright | 



I -Aright | 



(3.16) 



where we used the fact that |t2i| — |*22 1 near the reso- 
nance as will be seen below in Eq. (3.20). Hence the max- 
imum amplitude of the standing wave inside the cavity 
is 

|2 



M 



max | A rig hte l/Ci:r + A left e 



-\k\x I 



4|< 22 | 2 (3.17) 



iLin 



with its average half of Eq. (3.17). 

Now we relate the quantity \tii\ to the transfer ma- 
trix T^ v \ Because of the recursion relation (2.28), the 
resonance condition (3.5) reads 

(3.18) 



t22 



.iC/3 



tl2 2 +C-^ 3 t 22 



at ( = Q = £ r + ir;,., where we used the fact t\ 2 — t 2 \ on 
the basis of Eq. (2.30). We transform Eq. (3.18) into 



p 2r, r /3 



'■'I 



(3.19) 



For a very small \rj r \ 
sides in the form 

2 



-r] T , we can expand the both 



12 



0(r? r 




+ OfaO . (3.20) 



Using the relation (3.4), or |*i2 (^ r /3) | 2 + 1 
we can reduce Eq. (3.20) into 



*22 



2 fori 



1*22(6/3)1' 



(3.21) 



(Note again the fact that r\ x < 0. 11,12 ) The enhancement 
factor given in Eq. (3.17) is thereby 

6 



M 



(3.22) 



and the enhancement factor on average is 3/ |r? r |. This is 
indeed consistent with the numerical results in Fig. 7; the 
enhancement factor in Fig. 7 (a) is given by Table I (a) 
as M = 2.5 x 10 6 and that in Fig. 7 (b) is given by 
Table I (b) as M = 7.8 x 10 7 . 

Finally, we make a note for discussions in the final 
section §9. Notice that the field distribution in Fig. 7 is 
almost symmetric, being independent of the direction of 
the incident wave; although the incident wave, in fact, 
comes from the left in Fig. 7, the difference between the 
field strengths on the left end and the right end of the 
Cantor set is negligible. This is a common feature of the 
resonance. As we argued below Eq. (3.5), we have only 



outgoing waves at the resonance, and hence the reso- 
nant wave function itself is either symmetric or anti- 
symmetric, depending on the parity of the state. The 
field strengths on the left and the right are therefore 
equal. This is strictly correct only at the resonance, or on 
the complex plane. In the present case, however, the res- 
onance pole is so close to the real axis that the field dis- 
tribution is almost symmetric at the resonant frequency 
on the real axis, as we showed in Eqs. (3.16) and (3.20). 

We stress here that this happens for the three- 
dimensional resonance as well. For the Menger sponge 
embedded in three spatial dimensions, each resonant 
wave function corresponds to a irreducible representa- 
tion of the cubic point group Oh- Some of the resonant 
states have outgoing waves of the same amplitude in all 
six directions. At the resonant frequency close to such 
resonant states but on the real axis, the scattered wave 
has a strong amplitude even in the direction perpendic- 
ular to the incident wave vector. 

4. Renormalization 

In this section, we introduce a renormalization proce- 
dure of the transfer matrix for the Cantor set. We remind 
the readers of the definition of the transfer matrix T 2 in 
Eq. (2.21): 

T 2 (() = e _i ^" T »/ 2 e i ^ r " Tz e i ^" T! '/ 2 = e'^™^ cos0+a x sm</>) 

(4.1) 

The exponential operators of a y on the left-hand side rep- 
resents the effects of the boundaries at x = X and x = Y. 
We spread the effect of the reflection at the boundaries 
over the entire medium 2 in X < x < Y. In other 
words, we renormalized the boundary effect into a uni- 
form medium and hence had the single exponential oper- 
ator on the right-hand side, which takes the form of prop- 
agation in a uniform medium, similar to the form (2.25). 
The matrix a z in the exponent of the right-hand side 
represents the straight propagation, whereas the matrix 
a x represents the reflection inside the medium 2. 

We are now going to exploit the above idea and rep- 
resent the transfer matrix for the Cantor set of any gen- 
eration in the form of the transfer matrix for a uniform 
medium. In other words, we renormalize the fractal struc- 
ture into an effective uniform medium. 

We first prove that the transfer matrix for an 
arbitrary v is represented in the form 

with appropriate coefficients and . First, the sym- 
metry of the transfer matrix, Eq. (2.30), gives the sym- 
metry of its exponent as 

* logT( 1 )(C) = log^^C) = logT«(C). (4.3) 

Next, the exponents of the transfer matrices T\ and T 2 
are traceless. Therefore we have 

TrlogT (1) (C) = logdetT (1) (C) 



21ogdetT 2 ( | 



= TrlogTi 



2TrlogT 2 



0. (4.4) 
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Equations (4.3) and (4.4) indicate that the exponent of 
the matrix is symmetric and traceless. We can simi- 
larly show from the recursion formula (2.28) that the ex- 
ponents of all transfer matrices are symmetric and 
traceless. An arbitrary symmetric traceless 2x2 matrix 
is expressed by a linear combination of <r z and a x . Hence 
we have Eq. (4.2). 

Next, we prove that the coefficients and A^ are 
real for real (. We first show it for v = 0, or for T 2 defined 
by Eq. (2.21). In this case, we have 

s^=ncos(f> and A ^ = in sin 0. (4.5) 

Equations (2.22) and (2.23) indicate that s (0) and r (0) 
are both real (assuming that the dielectric constants are 
real). Next, we derive that s^ +1 ^ and A v+1 ^ are real from 
the assumption that and A^ are real. For the sake 
of the proof, we rotate in the recursion formula (2.28) 
the x and y axes of the spin space around the z axis by 
90 degrees and consider the identity 



Paper 

Equations (2.22) and (2.23) yield 



Author Name 



(»+D r 



r ("+l) 



<7y) 



<Ty)/3 e iC^/3 e iC( S ( " , ^-irC' ( T H )/3 _( 4 g) 



iC(( S ( " + 1) )*^-i(r(- + 1 ')V H ) 

= e -K( S '- v) a z -ir^a y )/3 e -Ka z /3 



xe 



<7y)/3 



(4.7) 



for ( e R. (Note here the fact that (a y )* = —<J y .) Hence 
we have 



e -ic((s < " +1) r^-i(r<"+ i >r ( T !/ ) = e 



-iC(s ( " + 1) ^-ir (, ' + 1) <T B ) 



(4.8) 

or, by taking the logarithm, 

(s^ 1 )) *a z - i{A u+ ^)*a y = s {v+ ^a z - w {v+ ^a y . (4.9) 

Because both sides of the equation are traceless 2x2 ma- 
trices, the expression with respect to the Pauli matrices 
must be unique. Therefore we have (s^ +1 ')* = s^ +1 ^ 
and (t> +1 ))* = A v+ ^ for ( £ R. 

To summarize this section, the transfer matrix given 
by Eq. (2.28) is expressed in the form 



TW(C) = e 1 ^ ^- ir ^ (4.10) 

with real coefficients and A v \ The former coefficient 
represents the (renormalized) straight propagation and 
the latter represents the (renormalized) reflection. 

5. Perturbation 

We now would like to find fixed points (s^*\A*^) of 
the renormalization transformation 



( s w,r w ) — ► (A- • *',r 



>+i) r (v+ih 



(5.1) 



As a first step, we here argue the transformation pertur- 
batively for e = e 2 je\ ~ 1. For simplicity, we hereafter 
assume /ii — /i 2 and hence n! — n. 
Our perturbation parameter is 

11=^- (5.2) 



A = 



1 fe 2 



s (0) = ncos</> = ^ ( — + 1 ) = 1 + A, (5.3) 



1 fe 2 



A ^ = in sin cp — 



1 =A. 



2 \si 

1 (e_2 

2 Ui 
Hence we have 

T(°)(C)^T 2 (C)=e i ^+ Ac(0) ) 

with 

A = o z and = a z — ia x . 



(5.4) 

(5.5) 
(5.6) 



Now we analyze the recursion relation (2.26), or the 
identity 

e iC(A+AC< >)/3 e iCA/3 e iC(A+AC<°>)/3 = e K(A+XC^ + 0(A 2 )) 

(5.7) 

using M. Suzuki's method of the "quantum anal- 
The calculation in Appendix A gives 



i) 13, 14 



ysis 
Eq. (A.18), or 



c« = ^-i J F( 1 )(CK, 



By taking the complex conjugate of Eq. (4.6), we have where 



F (D (C) 



X _ c iC/3 + e -iC/3 



Comparing the both sides of 



• C 2C • c 

sin i cos sin ^ 

3_ _ 2 3 3_ 



sin£ 



we have the renormalization transformation 
s (0) = 1 + A and A 0) = A — > 



(5.8) 

(5.9) 
(5.10) 

(5.11) 



s^-l + ^A and r (1) ~ AF^(C). (5.12) 
Similar calculation yields Eq. (A. 23), or 

s^~l+(J^j A and A 2) ~ AF^(C) (5.13) 



with 

F (2) (C) = F (l) (C)j p(l) 

We go forward similarly, arriving at 



(5.14) 



~ 1 + A and A V \Q AF^(C) (5.15) 



with 



f {v \c) = fW(c)^- 1 ) 



k=0 



= T 3 



2C 

sine ±L o K 
s k=i 



n 



We also define 



F(°)(C) = 1, 



(5.16) 



(5.17) 



which makes Eq. (5.15) hold for v = as well. 
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Fig. 8. Wc plot here (3/2)"FM (£) sin C for 1 < v < 5. Curves 
for v > 6 would overlap with the curve for v = 5 and be indis- 
tinguishable in this range. 



In the limit v — > oo, we have — ► 1 and the first- 
order perturbation vanishes, while the coefficient de- 
pends on £ oscillatorily. Numerical calculation shown in 
Fig. 8 suggests that the function F^(() has a limit of 
the form 

/<%)= lim (ff FW(C)sinC, (5.18) 

although we do not know its compact expression. For ( 
at which the function (£) / sin £ has a moderate value, 
the coefficient r^ v > decreases in the form 







(5.19) 



as v — > oo. The coefficient A u \ however, grows as 



r {v)(Q = 2 v \ for ( = ymn/2, 



(5.20) 



where m is a positive integer. This enhancement of the 
reflection amplitude comes from repeated reflections at 
the boundaries of the Cantor set. 

In order to highlight the cause of the enhance- 
ment (5.20), we here compare the above result for the 
Cantor set with the transfer matrix for an object in which 
the dielectric sheets are layered periodically. If the sys- 
tem of length L is divided into 21 + 1 alternating sheets 
of the medium 1 and the medium 2, the transfer matrix 
is given by 



C 



21 + 1 



Ti 



x T 2 



C 



21 + 1 

c 



21 + 1 



(5.21) 



where again £ = k\L. The Trotter formula gives its limit 
I — > oo in the form 



with 



y(°°) — c iC(so"z-ir(T x ) 
1 + S (°) 

and r = 



(5.22) 



— . (5.23) 



r (0) 

2 



Here we have no enhancement of the reflection amplitude 
such as Eq. (5.20). Indeed, the enhancement (5.20) oc- 
curs because the wave of ki = '& v miT /2L resonates in the 
cavity of the medium 1 of length L/3 as well as L/3 2 , 



L/3 3 and so on (Fig. 2) and the resonating waves inter- 
fere positively in every cavity. This never happens in the 
periodically layered medium (5.21). 

We stress here that the situation causing the strong 
resonance is common to the Menger sponge in three spa- 
tial dimensions; in both the Cantor set and the Monger 
sponge, the central part of the divided three segments 
is always a cavity. We thereby suggest that the recent 
experiment by Takeda et al. 1 found a strong resonance 
in the Menger sponge. 

6. Wave propagation within the perturbation 

We now compute the transmission and reflection co- 
efficients as well as seek the resonant states within the 
above first-order perturbation. Let us first expand the 
transfer matrix with respect to r^/s^ because 
sM = 0(1) and r-M = O(A) within the perturbation 
theory. Using again the quantum analysis described in 
Appendix A, we have Eq. (A. 35), or 



'12 



t22 



-Kb 



2 sM' 



sin 



-( s (")e- iCsW 



(6.1) 



(6.2) 



We thus arrive at the transmission coefficient T for ( e 
in the form (3.3) as 

1 



1 



,+^sin 2 ( Cs M) 



~ 1 



»(")■ 

7^r 



■sin 2 (0 M ) 



(6.3) 



Similarly we have the reflection coefficient R for ( e 

as 

i 2 



R = T x \Ti 



12 



sin 2 (C«M) 



(6.4) 



Within this order, we indeed have the flux conservation 
T + R = 1. With the use of Eq. (5.15), we have the 
expansion in A as 



R ~ 



~ A 



A 2 FM( C ) 2 xsin 2 
FM(C)sinC] 2 



c + 



= 2 2 ^A 2 sin 2 ^fTcos 2 ^. 

3«> 11 3*: 



k=l 



(6.5) 
(6.6) 



and T = 1 — R. Thus the reflection coefficient is enhanced 
as R ~ 2 2,y A 2 for £ = 3 l/ rmr/2, where m is a positive 
integer, whereas it vanishes for £ that gives the function 
/(*'(£)/ sin£ a moderate value. 

Incidentally, the sinusoidal function in the denomina- 
tor of F^(() in Eq. (6.5) is canceled by the factor sin£ 
and hence the function R does not diverge at £ = rmr. 
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Fig. 9. The reflection coefficient R due to the exact expres- 
sion (3.3) (the solid curves) and that due to the perturbational 
expression (6.7) (the gray curves) are plotted for < v < 4 
(from the top down to the bottom) and for (a) e = 1.01 and 
(b) e = 1.1. The difference between the solid and gray curves is 
actually indiscernible in the panel (a). Note that the horizontal 
axis of each panel is scaled by the factor 3". 



For the same cancellation to occur in Eq. (6.4), it may 
be better to modify the expression to 



R 



sin 



(6.7) 



so that the factor sin (C,s^ ) may cancel the denominator 
of the function (Cs^)- The difference between the 
expressions (6.4) and (6.7) lies only in higher orders of 
A, but we found that the latter expression (6.7) in fact 
reproduces the exact value of R better, particularly when 
the perturbation parameter A gets larger. We compared 
in Fig. 9 the values of R due to the exact expression (3.3) 
and due to the perturbational expression (6.7). The dif- 
ference is indiscernible for e = 1.01. A slight difference 
appears at length for e = 1.1 and v > 2. 



Fig. 9. Continued. 

We now look for the solution of the resonant condi- 
tion (3.5), or 



-its'-' 



2 S ^y 



sin 



) c -iC^> 



0. (6.c 



Rewriting the equation, we have 



-2iCs 



(") 



r W 2 i r w 2 



4 ' 



4sM 2 2sM 2 



(") 



4 ' 



(6.9) 



We thus arrive at the equation for the resonance pole £ r 
in the form 



2s(") 



(6.10) 



The use of Eq. (5.15) gives the expansion in A as 

A 



e-^ 2 Y M (&+^). 



(6.11) 
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(a) 




Fig. 10. (a) Three-dimensional plots of | (f ) | on the lower half 
plane of f for e = 1.01 and for < n < 4 (from the top down 
to the bottom), (b) Each plot of \F^(()\ is superimposed with 
the function e Imf = e - ^!. 



where Cr = 6- + i?7r- (Note, however, that the resonance 
poles exist in the lower half plane: r\ T < 0. 11,12 ) We im- 
mediately have 



and 



tan£ r = — 



ImfM^r + ir/r) 



(6.12) 



(6.13) 



RefM^r+if/r)' 

We can understand from Eq. (6.12) the fact that some 



Fig. 10. Continued. 



of the resonance poles approach the real axis for large 
v. For A = 0, the solution of Eq. (6.12) is r/ r = — oo. As 
we turn on the perturbation parameter A, the solution r] r 
takes a negatively large value. As we stressed above, how- 
ever, the function |-F^(C)| i s enhanced very much for 
some specific values of £ as is exemplified in Fig. 10 (a). 
For such values of £, the solution ?7 r of Eq. (6.12) ap- 
proaches the real axis as is shown in Fig. 10 (b). 

Next, we show that most of the resonance poles, upon 
climbing up to the real axis, sway either to the right or 
to the left as is schematically shown in Fig. 3. Because 
= 1, or ImF^ = 0, the resonance poles are aligned 
regularly at £ r = mir (m = 1, 2, 3, • • • ) for v = 0. For 



14 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



v = 1, the definition (5.9) gives 
F«(e + ir?) = 



e i£-V _ c i«/3-r)/3 + e -it/3+ V /3 _ e ~i£+V 



(6.14) 



e 1 ?-^ — e~^ +I ' 

which is, as r\ — > — oo, reduced to 

F«(£ + i77)~l-c- 2i ^ 3 + 2 "/ 3 

= (l - e 2 "/ 3 cos ^ + ie 2 "/ 3 sin ^, (6.15) 

and hence Eq. (6.13) is reduced to 

tan£ r ~ - c - 2 l*l/ 3 sin%. (6.16) 



The pole at £ r = tt for v = sways to the left for v = 1 
because Eq. (6.16) gives tan£ r < 0. Similarly, the pole at 
£ r = 2n for v = sways to the right for z/ = 1 because 
of tan£ r > 0, whereas the pole at £ r = 3ir does not sway. 
The same analysis gives 



-2i(|+i»,)/3 



p -2i(£+ir,)/9 



-2i(€+i»j)/9 



(6.17) 

= (l - e 2 "/ 9 cos + ie 2 "/ 9 sin f , (6.18) 



or 



tan£ r 



-2M/9 sin ^£_ 



for = 2 and gives for an arbitrary ^, 



tan£ r ~ -e- 2 l*l/ 3 %in 



26 
3^ 



(6.19) 



(6.20) 



The exponential factor in the right-hand side of 
Eq. (6.20) gets larger as v — ► oo. Remember that the 
amplitude of the function F^ v \x) has a peak centered 
at £ = 3"ir/2. Therefore the resonance poles to the left 
of the peak (where the poles climb up to the real axis) 
sways to the further left (tan 6 < for £ r < 3"7r/2), 
while the poles to the right of the peak sways to the fur- 
ther right (tan 6 > for £ r > 3 ly 7r/2). Thus we have the 
movement of the poles as is shown in Fig. 3. 

To summarize, the increase of the function (£) in 
some areas of £ causes the increase of the reflection co- 
efficient R. It also makes the resonance poles in the area 
approach the real axis, or makes the lifetimes of the res- 
onant states longer, and at the same time repulses the 
resonance poles away from its peak. Thus we end up with 
a wide area of a large reflection coefficient, beside which 
sharp peaks of the transmission coefficient appear. We 
here observed the above behavior in the framework of 
perturbation for e ~ 1, but we see the same behavior for 
e > 1 in Fig. 5 (e) and (f). We are going to argue below 
that the above behavior persists for e 3> 1 and for v 3> 1. 



7. The strongly renormalized limit 

In the previous section, we showed in the weakly renor- 
malized limit |s^| S> |r^| that the coefficient 
grows in the form (5.20) for some values of £ as the renor- 
malization proceeds, and accordingly the reflection coef- 
ficient R grows exponentially as is confirmed in Fig. 5 (a) . 
In the present section, we analyze the strongly renormal- 
ized limit |s^| <C and hence expand the transfer 
matrix with respect to s^/A u \ The computation de- 
scribed in Appendix A gives Eq. (A. 39), or 

<2 



12 



sinhCr M + 



2 7» 



■sinhCr M 



2 7» 



r Cr M coshCr M 



(7.1) 



To: 



cosh Cr {u) - ^sinhCr (,y) 



rCr M sinhCr M 



(7-2) 

2rM z 

The transmission coefficient T in the present limit is 
therefore given by 

T = 1 



T22*?22 



cosh 2 Cr M 



sinh 2 



r (i/) : 
1 



■Cr M sinhCr (,y) cos\\C,r {v) 



cosh 2 Cr(") 



1 + 



+ - 



■ tanh 



(r^ (Cr M -tanhCr (l/) ) 



(7.3) 



The reflection coefficient R is similarly given by 
R = Tx\T 12 \ 2 
~ Tsinh 2 Cr M 



+ ^ (l-t> (l/) cothCr (l/) ) 

r M 2 V s / 



cosh 2 Cr(") 
xtanhCr (,y) (( 



sinh 2 (r (u) - 



r" -tanhCr (l/) ) 



(7.4) 



Indeed we have the flux conservation T + R = 1 in this 
order of perturbation. For large |r^|, we have T — > 
and i? — > 1. 

The resonance condition T 2 2 = now reads 

>2 



cothCrM 



r(»0 



2rM ; 



(7.5) 
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By ignoring the second term of the right-hand side and 
rearranging further, we arrive at 



3 20 



(") 



1 + 



IS 



1 - 



is 



1 + 



2isW 



or 



3 (£+irj)r 



*i 1 + 



IS 



(") 



where we put £ 



£ + it]. We hence have 



1, 



or 



0. 



(7.6) 



(7.7) 



(7.8) 



which appears to be inconsistent with the first assump- 
tion | <C | . This means that the resonance poles 
do not exist in this limit. 

Thus we again observe the following behavior: in the 
area where grows, the poles move toward the real 
axis, but at the same time avoid the area of large \r^\, 
opening up a wide gap free of resonance poles. The re- 
flection coefficient R is close to unity in the wide gap, 
while the resonance poles just beside the gap give very 
sharp peaks of the transmission coefficient T. 

8. Trivial fixed point 

We showed in §5 that the reflection amplitude (£) 
grows exponentially for some £ as in Eq. (5.20), but de- 
cays exponentially for other £ as in Eq. (5.19). In the 
previous section, we analyze the limit of the growing re- 
flection amplitude. We examine in the present section, 
on the other hand, the trivial fixed point 



(0 



o. 



(8.1) 



We showed in §§6 and 7 that some of the resonance 
poles climb up to the real axis of the wave number near 
a peak of the reflection amplitude but are eventu- 
ally repulsed away from the peak. We argue here that 
such resonance poles in fact approaches the trivial fixed 
point (8.1). 

The transfer matrix becomes diagonal at the point 
where rM(£) = 0: 



-iCs 1 



which is immediately followed by 

T = 1 and R = Q. 



(8.2) 



(8.3) 



Thus the transmission coefficient has a sharp peak at the 
real solutions of the equation r^)(£) = as is demon- 
strated in Fig. 11. We also know that transmission co- 
efficient has a sharp peak in the case that a resonance 
pole is close to the real axis. Combining the above two 
observations, we conclude that the resonance pole, 

(1) first climbing up a peak of toward the real axis, 

(2) next sliding off the peak sideways, 

(3) eventually approaches a solution of the equation 
j-W (£) = on the real axis. 




Fig. 11. The bottom graph of each of the panel (a) v = 1 and 
the panel (b) v = 2 gives the f dependence of (the broken 
curves) and (the solid curves) for e = 10. The mid graph of 
each panel shows the locations of the resonance poles and the top 
graph of each panel shows the £ dependence of the transmission 
and reflection coefficients. The mid and the top graphs are the 
same as the ones in Fig. 5 (e) and (f). The vertical gray lines 
indicate the zeros of rW. Some of the zeros arc common to the 
both panels. 



We indeed observe in Fig. 11 that the resonance poles 
right next to the wide regions of R <~ 1 are particularly 
close to the real axis and are accompanied by a very 
narrow peak of T . The poles away from the regions R ~ 
1, on the other hand, arc not close to peaks of A u \ and 
hence are not close to the fixed point — on the real 
axis. (We have T — > oo at the resonance pole, while we 
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have T = 1 on the real axis, which appears to lead to 
a contradiction when the resonance pole approaches the 
real axis. They are reconciled if ds^ v \C) / dC, diverges at 
the pole.) 

Incidentally, the recursion formula of the transfer ma- 
trix shows that the matrix T^ V+1 \Q) is diagonal if the 
matrix TM(C/3) is diagonal. This means r^ +1 )(C) = 
if r^(C/3) = 0. This is also demonstrated in Fig. 11; 
some of the solutions r^(() = are indeed given by 
r«(C/3) = 0. 

9. Discussions 

To summarize, we have analyzed how the resonance 
poles Cr move in the complex wave-number plane as we 
progress to higher generations of the Cantor set. We con- 
clude that the fractal structure of the medium produces 
resonance poles very close to the real axis of the wave 
number with a wide region of the wave number without 
resonant states. As a result, we have sharp peaks of the 
transmission coefficient (T <~ 1) and a wide region of 
nearly total reflection (R ~ 1) between them. 

Although we carried out our analysis only in one spa- 
tial dimension, we expect that the above conclusion itself 
is independent of the dimensionality. As we emphasized 
in Fig. 2 and below Eq. (5.20), the strong resonance is 
essentially due to the fractal structure in which the cen- 
tral part of the divided three segments is always a cavity. 
This structure is common to the Monger sponge in three 
spatial dimensions. We thereby suggest that the recent 
experiment by Takeda et al. 1 found a strong resonance 
in the Menger sponge. 

It is then quite natural for them to observe dips both 
in the transmission amplitude and the reflection ampli- 
tude. In the resonant state, the incident wave is scattered 
strongly and repeatedly inside the cavity for the duration 
of the resonant lifetime. Therefore the outgoing wave can 
have an amplitude over quite a wide range of the scatter- 
ing angle in three dimensions. We remind the readers of 
the fact that, in one spatial dimension, the resonant wave 
function itself is either symmetric or anti-symmetric and 
the field strengths on the left and the right are therefore 
equal. Close to the resonance pole, the field distribution 
is almost symmetric at the resonant frequency on the real 
axis. Indeed, the resonant scattering observed in Fig. 7 
is almost symmetric and independent of the direction of 
the incident wave. We stress here that this happens for 
the three-dimensional resonance as well. For the Menger 
sponge embedded in three spatial dimensions, some of 
the resonant states have outgoing waves of the same am- 
plitude in all six directions. At the resonant frequency 
close to such resonant states but on the real axis, the 
scattered wave has a strong amplitude even in the di- 
rections perpendicular to the incident wave vector. We 
hence claim that the lost fraction of the transmission and 
reflection amplitudes in the experiment by Takeda et al. 1 
must have been emitted in the directions normal to the 
incident direction. 
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Appendix A. Suzuki's quantum analysis 

In the present Appendix, we give details of the calcula- 
tions producing Eqs. (5.8), (5.13), (6.1)-(6.2) and (7.1)- 
(7.2). We heavily used M. Suzuki's "quantum analy- 
sis" 13, 14 in the calculation. See Hatano and Suzuki 15 for 
a recent review. 

The main feature of the quantum analysis is the differ- 
entiation of an operator functional f(A) by an operator 
A. The derivative df{A)/dA is given as a hyperoperator 
mapping dA to df(A), defined in 

df(A) = hm W + Mf-M s M . ^ (A.1) 
v ; h^o h dA K ' 

where the convergence of the limit operation is in the 
sense of the norm convergence uniform with respect to 
the arbitrary operator dA. It should be noted that the 
hyperoperator df(A)/dA is expressed in terms of the 
operator A itself and its inner derivation 6a- The in- 
ner derivation is a hyperoperator giving the operation of 
nothing other than the commutation relation: 

6 A B=[A,B}. (A.2) 

Once we obtain the operator derivative (A.l), we can 
calculate the parameter derivative of the operator func- 
tional f(A(x)) in the form 

df(A(x)) _ df(A) dA(x) 



dx 



dA 



dx 



(A.3) 
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A particularly important formula given by the quantum 
analysis is the operator derivative of the exponential of 
an operator: 14 
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dc A 

~dA 



1 



^—Sa 



= e 



5a 



71=1 



5a 



n-l 



(A.4) 



We use this formula repeatedly hereafter. 

Let us differentiate the both sides of Eq. (5.7) with 
respect to A. The formula (A.4) with Eq. (A. 3) gives 

d iC04+AC<°>)/3 



dA 



= [ A p K(A+\C^)/3 



1 - cxp [S iaA+xcm)/3 ] g(0) 



J iC(A+AC<°))/3 

Note that 5 x a = x5a- Thus we have 



5) 



dA 



(iC/3)<5A 



A=0 



(K/Wa 



C (0) 



5a 



(A.6) 



XA/3t ± C (0) e -iCA/3 e i(A/3 



Sa 

= c (iC/3U, e-^-l 



5a 



C (0) e KA/3 



1 _ C (K/3)S A 

T A 



C r(0) e iCA/3_ 



Here we used the identity 

e^Be"^ = e^B 



(A.7) 



(A.8) 



in moving from the third line to the fourth line. There- 
fore, the differentiation of the left-hand side of Eq. (5.7) 
yields 

d 



dA 



e iC(A+AC< >)/3 e iCA/3 e iC(A+AC< >)/3 



A=0 



1 _ e (K/3)S A 



_ e KA 



5a 

e ~(K/3)S A _ 1 

5~ A 



C(o) e i^ 

C (0) , (A.9) 



where we used both of the expressions (A.6) and (A.7). 
On the other hand, the differentiation of the right-hand 
side of Eq. (5.7) yields 

-KS A 



^_AC(A+\C {1) +0(\ 2 )) 

dA 



-C (1) (A.10) 



A=0 S A 

because of the formula (A.4). Comparison of Eq. (A.9) 



5a 



-i(6 A 



-i(A 



- 1 

,(iC/3)«5/ 



5 A 



C (0) e i<A 

e -(K/3)5 A _ 1 



+e 



i(Al 



-KS A 



- 1 
+ 



C (0) 

5a 

(i _ C (K/3)S A ^ c (o) 



C -KS A _ e -(2iC/3)5 A + e -(K/3)S A _ 1 



C (0) , (A.ll) 



g-iC<5A _ 1 

where we again used Eq. (A.8) in moving from the first 
line to the second line. The final expression is an even 
function of (5a- Therefore the Taylor expansion with re- 
spect to (5a gives only even-order terms: 

C (1) = ]T c k (-K) k 5 A k C (a \ (A.12) 

/c:cvcn 

where the coefficients c k are given by the Taylor expan- 
sion 

_ p 2x/3 + e x/3 ■ 



CkX 



(A.13) 



with Co = 2/3. 

Using the definition (5.6) explicitly here, we have 

S A C^ = 



[a z , a z — ia x \ 

= — i [o'zjO'x] = 2 a y, 

5 A 2 C< 0) = [a z ,2a y ] = -4i<r x , 
5 A 3 C^ = [a z ,-4ia x ] = 8a y ,- 



or 



6 A k C^ = ! f ay 

{ O ka <T z 



for odd k, 
2 k \a T for even k. 



(A.14) 
(A.15) 
(A.16) 

(A.17) 



Substituting the result (A.17) for the terms in the sum- 
mation of Eq. (A.12), we arrive at 

\<J Z - i<?x ^2 c fe (-2iC) fc 



C (i) = 



fcicvcn 



2 e" 2 * - e- 4i ?/ 3 + e- 2i ?/ 3 - 1 

r z ^r—i 1<Jx 



= -v,-iFW(Qv x , 



(A.18) 



or Eq. (5.8), where F (1) (C) is defined by Eq. (5.10). 

The next renormalization transformation is obtained 
by replacing with and C« with C< 2 ) in 

Eq. (A.ll). In the replaced the argument of the 

function is C/3. Instead of Eqs. (A.14)-(A.17), we 
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(A.19) 



6 A 2 C^ 



S A 3 C^ = 



-4^(0., 



(A.20) 



(A.21) 



or 



5 A k CW 



2 
3' 



for odd fc, 

|4o^ - 2 fe LF (1) ^0 cr* for even k. 

(A.22) 

We thus have 

= (0 2 ^-^ (i) (o^ (i) (0^ 

= (0 ^-iF( 2 )(CK, (A.23) 
or Eq. (5.13), where F( 2 )(C) is defined by (5.14). 
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We now expand the transfer matrix with respect to / s^ v \ Calculation similar to Eqs. (A. 7) yields 



d :^„^) ax) 1 ~ exp [-iC i/ v K. - irM5 a .)] 



drW 



(&x 



fc=l 



k-1 



(A.24) 



The zeroth-ordcr term of in the summation on the right-hand side is 

E \-KS y ') x fe_l 
n d <?z a x- 



fc=l 



k\ 



The same computation as in Eqs. (A.14)-(A.16) gives 

2 fe icr„ 



5<T, Ox 



2 k a x 



for odd k, 
for even k, 



(A.25) 



(A.26) 



and hence we have 



~ (-2iCs^) fe - 1 A (-2iC S M) fc - 1 

Eq. (A.25)= Yl \, ^+ E - 



fc=l,3,5,- 



fc=2,4,6,- 



fc! 



2Cs( 



^- |a 2; sin2Cs (l/) +cr y l-cos2Cs M } 



The first-order term of in the summation on the right-hand side of Eq. (A.24) is 



i £ - M r<"> £ ^ fe "'" 2 ^ • 



fe-i 



fc-2 



fc=2 



fc! 



Z=0 



Because of Eq. (A.26), we have 







for odd Z, 
for even Z. 



1 -cos2Cs M 

(A.27) 

(A.28) 
(A.29) 



Therefore, the only remaining terms in the summations of Eq. (A.28) are for I = k — 2 = odd with the result 
5 a J az k - 2 a x = -2 k - x o z . Thus we have 



Eq. (A.28) = ir» £ 

fc=3,5,7,- 
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To summarize, we have 



and then have 
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(A.30) 



(A.31) 
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We finally obtain 
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(A.34) 
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(A.35) 



(A.36) 



which gives Eqs. (6.1) and (6.2). 

Finally, we expand the transfer matrix with respect to s^/r^ in the limit opposite to the above. We can actually 
use the expression (A.36) of the expansion with respect to /s^K Rotating the axes of the spin space around the 
spin y axis, we have the transformation 

<r z — >o x , a x — >-a z , and hence T H = e K(s^a z -ir^a x ) — > f( v ) = jW^.+ir™*,) ( A3? ) 

We can thereby obtain the expansion with respect to s^/A^ by replacing and in Eq. (A.36) with ir^ 
and is^\ respectively. The result is 



j-W ~ e -Cr 



<Jz- 



(A.38) 



Then we rotate the axis back in the way opposite to Eq. (A. 37), arriving at 
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which gives Eqs. (7.1) and (7.2). 
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